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ABSTRACT 

In this second paper of the series, we present the results from optical Gemini-North GMOS-IFU and 
WIYN DcnscPak IFU spectroscopic observations of the starburst and inner wind zones of M82, with 
a focus on the state of the T ~ 10 4 K ionized interstellar medium. Our electron density maps show 
peaks of a few 1000 cm' 3 (implying very high thermal pressures), local small spatial-scale variations, 
and a fall-off in the minor axis direction. We discuss the implications of these results with regards to 
the conditions/locations that may favour the escape of individual cluster winds that ultimately power 
the large-scale superwind. 

Our findings, when combined with the body of literature built up over the last decade on the state 
of the ISM in M82, imply that the starburst environment is highly fragmented into a range of clouds 
from small/dense clumps with low filling factors (<1 pc, n e >10 4 cm -3 ) to larger filling factor, less 
dense gas. The most compact clouds seem to be found in the cores of the star cluster complexes, 
whereas the cloud sizes in the inter-complex region are larger. These dense clouds are bathed with an 
intense radiation field and embedded in an extensive high temperature (T>10 6 K), X-ray-emitting 
ISM that is a product of the high star formation rates in the starburst zones of M82. The near-constant 
state of the ionization state of the ~T0 4 K gas throughout the M82 starburst zone can be explained 
as a consequence of the small cloud sizes, which allow the gas conditions to respond quickly to any 
changes. In Paper I we found that the observed emission lines are composed of multiple components, 
including a broad (FWHM ~ 150-350 kms" 1 ) feature that we associate with emission from turbulent 
mixing layers (TMLs) on the surfaces of the gas clouds, resulting from the interaction of the fast wind 
outflows from the SSCs. The large number of compact clouds and wind sources provides an ideal 
environment for broad line emission, and explains the large observed broad/narrow line flux ratios. 

We have examined in more detail the discrete outflow channel identified within the inner wind in 
Paper I. The channel appears as a coherent, expanding cylindrical structure of length >120 pc and 
and width 35-50 pc. The walls maintain an approximately constant (but subsonic) expansion velocity 
of ~60 kms -1 , and are defined by peaks and troughs in the densities of the different line components. 
We hypothesise that as the hot wind fluid flows down the channel cavity, it interacts with the cooler, 
denser walls of the channel and with entrained material within the flow to produce broad-line emission, 
while the walls themselves emit primarily the narrow lines. We use the channel to examine further the 
relationship between the narrow and broad component emitting gas within the inner wind. Within the 
starburst energy injection zone, we find that turbulent motions (as traced by the broad component) 
appear to play an increasing role with height. 

Finally, we have argued that a point-like knot identified in GMOS position 4, exhibiting blueshifted 
(by ~140 kms -1 ), broad (<350 kms -1 ) Ha emission and enhanced [Sn]/Ha and [Nil] /Ha ratios, is 
most likely an ejected luminous blue variable (LBV) type object. 

Subject headings: galaxies: starburst - galaxies: evolution - galaxies: individual: M82 - galaxies: ISM 



1. INTRODUCTION 

Understanding feedback effects between massive stars, 
star clusters and the interstellar medium (ISM) is fun- 
damental to the study of galaxy evolution. In this sec- 
ond paper of the series presenting high resolution opti- 
cal integral field unit (IFU) observations of ionized gas 
in the starburst galaxy M82, we focus on the state of 
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the nebular gas found in the disk and inner wind of 
this galaxy. This follows on from a compan ion paper 
(| Westmoquette et al.ll2009l hereafter iPaper If ) where we 
presented the data, described our observation and data 
reduction methods, and discussed the gaseous and stel- 
lar dynamics of the starburst core and inner disk of the 
galaxy. 

M82 has long been known to be one of the nearest 
starburst galaxies with a large-scale optica lly-bright su- 
perwind (see, e.g.. lLvnds fc Sand agc 1961). This makes 
it of fundamental importance as an analogue to in- 
tensely star- forming galaxies at high-z. The current 
(~10 Myr) starburst activity is centred on the nucleus 
with a diameter of ^500 pc (~30"). Within this re- 
gion there are a number of prominent, high surface- 



2 



M. S. Westmoquette et al. 



brightness complexes 7 (labelled A-E), each known to 
contain many hu ndre ds of young massive star clusters 
(jO'Connell et al.lll995l : lMelo et al.ll2005t) . It is the energy 
from these clusters that drives the outflow to kpc-scales 
(IShopbell fe Bland-Hawthoral[l998t IStevens et aD [20031 : 
Stricklan d fe Heck man 2003). 

In iPaper II we found that the ionized gas within the 
starburst is dynamically very complex with much lo- 
calised line splitting, and many overlapping expanding 
structures located at different radii. The high signal- 
to-noise of the data were sufficient for us to be able to 
accurately decompose the emission line profiles into mul- 
tiple narrow components (FWHM ~ 30-130 kms -1 ) su- 
perimposed on a broad (FWHM ~ 150-350 kms -1 ) fea- 
ture. F ollowing the recent work by IWestmoquette et al.l 
(l2007al lbllcl. the latter will hereafter be referred to as 
lW07aV we identified the broad component as originat- 
ing from turbulent layers (|Slavin et al.lll993h on the sur- 
faces of interstellar clouds set up by the impact of high 
energy photons and fast-fl owing winds from th e embed- 
ded young star clu sters (|Pittard et al.l 120051 . see also 
iBinette et al.l [2009) . These layers may also be where 
mass is also stripp ed and entrained an d/or loaded into 
the win d outflow (|Suchkov et al.|[i~994}) . Results from 
IPaper II show that turbulent gas mixing into outflowing 
hot gas becomes increasingly important with height be- 
low the main disk of M82. Furthermore, careful inspec- 
tion of the ionized gas properties directly below complex 
C revealed a discrete channel within the outflow. Fur- 
ther investigation of this interesting structure forms part 
of the focus of this paper. 

High central densities, n , (hence thermal pressures) 
and steep radial density gradients have l ong been known 
in M82 (lO'Connell fe Mangand 119781 : iHeckman et~all 
119901: iMcKeith et al.lll995l ), and formed one of the key 
pieces of evidence supporting the starburst/superwind 
explanation of the M82 "explosion" (first reported by 
iLvnds fc Sandagelll963f ). On the most basic level, the 
collective effect of the supernovae and stellar winds in a 
starburst blows a bubble of hot gas that is driven out- 
wards by a pressure much greater than the surrounding 
ambient medium. The overpressurisation of this bubble 
drives it to expand further and, if the energy input is suf- 
ficient, to blow-out of the galactic disk, allowing a coher- 
ent outflow of material to develop. The absolute value of 
the central density measured in M82 has, to some degree, 
been a function of the spatial resolution of the observa- 
tions, since the nuclear regions of M82 exhibit a complex 
morphology of complexes, voids and dust lanes. From 
HS T /STIS spectroscopy, w e meas ured peak densities in 
the range 1500-2000 cm" 3 (|W07aY Until now, densities 
have mostly been measured using one-dimensional slits 
(albeit at various positions and PAs), but this method 
is clearly limited in the amount of information available 
regarding the full two-dimensional variation of the den- 
sities. 

Spatially resolved diagnostics of the gas excita- 
tion state, however, are easier to come by since, 
in the crudest form, two narrow-band images can 

7 In this w ork we adopt the term starburst complexes 
IjBastian et al.l [2005) to describe regions of concentrated 
star/cluster formation in the M82 disk. We prefer this term to 
the one of starburst clumps used in Paper I since this avoids con- 
fusion with gas clumps. 



simply be ratioed. 
measured from 

throug h to far-IR at a number of 
tions (ISatyapal et al.l 119951: lAchtermann fc Lacyl 199 
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a number of spatial resolu- 
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The dominant finding of these studies is that there is very 
little variation in the line ratios within the starburst re- 
gion, indicating a high degree of uniformity of conditions 
within the ionized gas. This is surprising considering the 
dramatic density variations and the existence of multiple, 
spatially defined, luminous starburst complexes. 

Formed only with the last decade, the emerging pic- 
ture of conditions within the M82 starburst has re- 
sulted from the comparison of multi-wavelength molec- 
ular and atomic gas observations with detailed PDR 
and photoionization codes. To summarise, the ISM is 
thought to be highly fragmented into a continuum of 
cloud sizes and densities (ranging from <lp c and <few 
x 10 4 cm~ 3 ), with l ow are a covering factors (lLord et al.l 
19961: IColbert et al.lll999l: iForster Schreiber et al.l 12001k 



20031 iRodriguez-Rico et al.ll2004t IW07cf ). A recent study 
bv lFuente et all (|2008l ) indicated that ~90% of the mass 
of the molecular gas exists within small clouds with 
Av > 5 mag. These clouds are highly pressured , 
and partly disru pted and dissociated (jMao et al.l 120001 : 
iWeifi et al.1 12001 as a result of being in such a high en- 
ergy environment. The clouds are embedded in a hot, 
diffuse surrounding mediu m of equal pressure (recently 
detec ted and modelled by IStrickland fc Heckmanl 120071 . 
l2009t ) maintained by the large number of SN shocks prop- 
agating through the starburst. In order to m aintain the 
observed levels of ionization and excitation, lLord et alJ 
found the average separation of the clouds and ioniz- 
ing sources had to be of order their size (1-7 pc). This 
agrees w ell with the avera ge cluster separations cata- 
logued bv lMelo et~aT1 (|2005h with HST observations, and 
with high-resolu tion CO observations (|Mao et al.l 120001 : 
IWeifi et aT1l200lD . Thus, the fragmented and well-mixed 
nature of the ISM within the starburst must allow it to 
rapidly equalise after environmental changes, ex plainin g 
t he aforementioned consistency of the line r atios (|W07cf). 



IShopbell fc Bland-Hawthornl (|1998f l and IBeirao et all 
(|2008| ) found evidence for changes in the gas excitation 
out of the plane along outflow axis (minor axis). This 
has been interpreted as resulting from a combination of 
a decrease in gas density leading to an increase in the 
ionization parameter, dilution of radiation field, and the 
increasing importance of shocks in the halo. Studying the 
relationship between the gas conditions in the starburst 
and the outflow, and the transition between the two, is 
of essential importance for understanding how the wind 
develops and evolves. 

In this paper we will focus on the derived nebular prop- 
erties of the ionized gas, including the densities and ex- 
citations. By combining both detailed (~0'.'7) observa- 
tions of the starburst core and inner wind made with the 
Gemini GMOS-IFU, and wider-field observations of the 
inner ~2 kpc of the disk made with the DensePak in- 
strument, we will attempt to build up a coherent model 
of the ISM conditions (physical composition, densities, 
excitations) within the starburst, use the density distri- 
bution to discuss implications for the escape of individual 
cluster winds in different parts of the starburst, and ex- 
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amine the ionized gas flow properties within the inner 
wind and the starburst energy injection zone. 

In this work, we have adopted a distance to M82 o f 
3.6 Mpc (jFreedman et alJll994tlMcCommas et al.ll2009l ). 
meaning 1" = 17.5 pc. 

2. OBSERVATIONS AND DATA REDUCTION 

A full description of our obser vations and data reduc- 
tion methods is given in lPaperll . In the following we will 
briefly recap some of the main points. 

2.1. Gemini observations 

In February 2006 and February 2007 observations using 
the Gemini-North Multi-O bject Spectrograph (CMOS) 
Integral Field Unit (IFU; lAllington-Smith et all l2002fi 
were obtained covering five regions near the centre of 
M82 (programme IDs: GN-2006A-Q-38, and GN-2007A- 
Q-21, PI: L.J. Smith), with 0.3-0.8 arcsec seeing (5-14 
pc at the distance of M82; see indicator on Fig. [2]). The 
IFU in two-slit mode gives a field-of-view of 7 x 5 arc- 
sees (which corresponds to approximately 50 x 35 pc at 
the distance of M82) sampled by 1000 hexagonal con- 
tiguous fibres of 0'.'2 diameter. A separate block of sky 
fibres, offset by 1 arcmin from the object field, provides 
simultaneous sky observations. We took two dithered ex- 
posures at each position with integration times between 
1200 and 2400 s each, and used the R831 grating to give 
a spectral coverage of 6100-6790 A with a dispersion of 
0.34 A pix . A spectral resolution of R ~ 5200 was 
measured from single Gaussian fits to isolated arc lines. 
The locations of each field are shown on an HST / ACS 
F658N image in Fig. Q] 

Basic data reduction was performed following the stan- 
dard Gemini reduction pipeline (implemented in iraf 8 ). 
After reduction and combining the dithered exposures 
the five final data cubes each contained 825 wavelength 
and flux calibrated, sky subtracted spectra. 

2.2. WIYN observations 

DensePak (|Barden et alj 119981 ) was a small fibre-fed 
integral field array attached to the Nasmyth focus of the 
WIYN 9 3.5-m telescope. It had 91 3" fibres arranged in 
an array of 30 x 45 arcsecs. Four additional fibres, offset 
by ^60" from the array centre, served as dedicated sky 
fibres. 

On 14th April 2001, we used DensePak to observe the 
inner disk of M82. The 860 line mm -1 grating was used 
in two settings; four fields were observed in the wave- 
length range of 5820-6755 A (dispersion 0.46 A pix -1 ; 
R ~ 6700), three of these positions were also observed 
in the wavelength range of 7745-9700 A (dispersion 
0.96 A pix 1 ; R ~ 3650). These settings allowed us to 
cover many of the brightest optical and near-IR (NIR) 
nebular emission and stellar absorption lines. A number 
of bias frames, flat-fields and arc calibration exposures 
were also taken together with the science frames. Basic 

8 The Image Reduction and Analysis Facility (iraf) is dis- 
tributed by the National Optical Astronomy Observatories which is 
operated by the Association of Universities for Research in Astron- 
omy, Inc. under cooperative agreement with the National Science 
Foundation. 

9 The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatories. 



reduction was achieved using standard IRAF packages. 
After reduction, the datafiles contained 86 wavelength 
calibrated and sky-subtracted spectra. Fig. [1] shows the 
position of the DensePak fields overlaid on the HST 
F658N image. 

2.3. Decomposing the line profiles 

Following the metho dology employed by 
IWestmoquette et all (|2007al lbl). we fitted multiple 
Gaussian profile models to each emission line u sing an 
idl Y 2 -fitting package called pan (Peak ANalvsis: iDimed 
l2005i) , to quantify the gas properties observed in each 
IFU field. 

Ha is detected in every fibre of both the GMOS and 
DensePak datasets. The high S/N and spectral resolu- 
tion of our data allowed us to quantify the line profile 
shapes to a high degree of accuracy. In general, we find 
the emission lines to be composed of a bright, narrow 
component (hereafter referred to as CI; this was some- 
times split into two narrow components, hereafter CI 
and C3) overlaid on a fainter, broad component (here- 
after C2). In some regions of our DensePak data, we can 
also identify an Ha absorption component of stellar ori- 
gin (hereafter C4). A full description of the procedure 
we employed to fit the lines and automatically deter- 
mine th e number of Gaussian models to fit is given in 
iPaper J . together with estimates of the fit uncertainties. 
It is prudent to recall that for the purposes of consis- 
tency in the fitting and analysis, where a triple-Gaussian 
fit was chosen we specified that the broadest component 
should be assigned to component 2 (C2), and after that, 
the brightest to be component 1 (CI), and the fainter 
narrow component to component 3 (C3). 

Given the limitations of current IFU observing tech- 
niques and multi-component line fitting methods, a rig- 
orous analysis of the uncertainties associated with the 
results presented in this paper is currently not feasible. 
Firstly, the flux calibration of IFU data is not just sub- 
ject to the usual kinds of uncertainties found in spectro- 
scopic data, but also to additional complications due to 
the complexity and design of the instrument. Further- 
more, the combination and interpolation of the dithered 
IFU maps introduces additional uncertainties. These 
factors together mean that we cannot rigorously quan- 
tify the uncertainties in the calibration of the individual 
spectra. Secondly, although a reduced % 2 value for a 
multi-component line fit can be calculated, we currently 
do not fully understand how to quantify the uncertainties 
associated with each individual Gaussian component fit 
to the line profile. However, we can estimate errors from 
the observed variances in observed parameters in regions 
where conditions a re relatiy ely stable. 

As discussed in IPaper fl . we estimate that, for the 
fluxes, there is a 0.5-10 per cent error in CI (with the 
range resulting mostly from the variation in S /N levels) , 
8-15 per cent on C2 and 10-80 per cent on the C3 flux, 
with the level depending on the relative flux in C3. The 
presence of a third fit component may slightly increase 
the errors in the CI and C2 fluxes particularly where C3 
is faint. When these estimates are carried through, we 
might expect similar percentage errors for the line ratio 
and electron density measurements. Due to the afore- 
mentioned difficulties, we do not attempt to include error 
maps or show error bars for each figure, but where spe- 
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cific measurements are quoted in the text we have made 
an effort to indicate approximate uncertainties. Finally, 
we have converted the fluxes of the Ha absorption com- 
ponent (C4) to equivalent widths (EWs), and we esti- 
mate that the uncertainties on the EWs lie in the range 
of ~ 0.5-l A 

In iPaper 1 we presented the line width and radial ve- 
locity results from the GMOS and DensePak datasets. In 
the sections following this, we will present and describe 
the remainder of our results, including the line flux and 
flux-derived properties of the emission line gas. 

3. NEBULAR PROPERTIES OF THE NUCLEAR REGIONS 
AND INNER WIND 

3.1. Ha flux distribution 

Figs. [2j [3] and [4] show maps of the Ha fluxes for the 
three line components, respectively. For reference, we 
have overplotted contours from the EST/ ACS F439W 
(-B-band; iMutchler et"aTll2007f) image in order to high- 
light the distribution of stellar emission and therefore the 
locations and extents of regions A and C. We have also 
indicated the maximum size of the seeing disk experi- 
enced during these observations (C/8) with filled circle. 
We begin by noting that CI and C2 are detected over 
almost all of the five positions, and C3 is identified in 
approximately one third of all spaxels. We also note 
that the typical EW of the Ha emission line in position 
1 is 50 A ' meaning that we can safely ignore the effects 
of underlying stellar absorption. 

To a first-order approximation the flux morphology in 
the three Ha line components is fairly similar. However 
there are important differences. Firstly, we notice that 
the region aro und the cluster M 82-A1 (northern corner 
of position 1; ISmith et al.l 12006^ 1 does not stand out in 
any of the three component maps (although it is clearly 
seen in the F439W contours). Secondly, the C2 and C3 
fluxes appear to be brighter on the southern side of any 
prominent feature compared to CI. For instance in po- 
sition 1 (region A) the flux of CI peaks to the north and 
east of region A, whereas C2 peaks in the centre to south- 
west of region A (its flux, however, is strongly influenced 
by the varying contribution from C3). Furthermore, the 
peak in CI flux coincident with the large Ha knot in po- 
sition 3 is located slightly to the north compared to the 
C2 peak. In position 2, this offset trend is again seen, 
with the peaks in C2 and C3 being offset slightly to the 
south compared to CI. 

In order to assess the importance of C2 in terms of its 
contribution to the total Ha flux, we have mapped the 
the C2/C1 flux ratio for all five IFU positions in Fig. [5l 
Surprisingly, of the spaxels that have both a CI and C2 
detection, ~60% exhibit C2/C1 flux ratios greater than 
unity (i.e. on average C2 contributes more to the total 
line flux than CI; see also Table [1]). Ratios >10 most 
likely result simply from bad fits. However, true ratios 
in the ran ge 5-10 are seen in the wind positions 4 and 5 
(fig. 4f of iPaperll shows an example line profile with a 
C2/C1 flux ratio in this range), and there are a significant 
number of spaxels with C2/C1 ratios of ^2 in position 3 
and parts of the south-west of position 1. The majority 
of spaxels with ratios <1 are found in positions 1 and 2 
(regions A and C), whereas those exhibiting the lowest 
ratios (0.1-0.3) are found in position 4 (coincident with 
the regions containing the broadest C2 emission, i.e. the 



position 4 knot - see lPaper J ). In summary, although the 
C2/C1 flux ratio varies significantly and rapidly over the 
five IFU fields, the general trend is of increasing C2/C1 
with radius. This is discussed further and in context with 
the rest of our findings in Section 15.31 Finally, we note 
that finding such a large number of high (>1) C2/C1 flux 
ratios i s in sharp contrast with wha t was found in NGC 
1569 bv lWestmoauette et al.l (|2007allbT ). There, ratios fell 
in the range 0.1-0.3 with an average 0.12. 

3.2. Electron densities 

Figs. [6j [7] and [8] show maps of the spatial variation 
of electron density, n e , in each of the three line compo- 
nents, derived from flux ratios between Gaussian fits to 
the [S n] AA6717,6731 lines (assuming an electron temper- 
ature, T c = 10 4 K). 

The density distributions in the three line components 
are distinctly different. We begin by discussing the densi- 
ties in CI (Fig. [6]). Within our limited spatial coverage, 
there are two peaks in n c , one in the northern corner 
of position 1 and one in the northern corner of position 
2. Position 1 covers the eastern half of region A; the 
peak in n c at ~2000-3000 cm~ 3 is coincident with the 
region around cluster M82-A1 and the well-known major 
axis dust lane that bis ects the whole nuclear region (see 
lO'Connell et al.lll995f ). From here, n e falls off rapidly 
to the south and west, meaning that in the core and 
western half of region A (position 3) the densities are 
significantly lower. The peak in position 2 is again lo- 
cated to the north of the complex core (region C), and 
reaches ^2000-3000 cm -3 . Following the same trend, 
the densities then fall off towards the south, but remain 
fairly constant at ^1000 cm -3 throughout the rest of re- 
gion C. The se values are consistent with previous studies 
(e.g. IW07d ). The densities in positions 3, 4 and 5 fall 
off consistently towards the south-cast (along the mi- 
nor axis direction) , almost reaching the low density limit 
(~100 cm -3 ) at the edge of our coverage. Interestingly, 
the density morphology in both the wind positions (4 and 
5) follows similar linear minor axis-oriented patterns as 
observ ed previ ously in the FWHM and radial velocity 
maps (IPaper J ). In position 5, we associated these with 
a discrete outflow channel, but finding these type of pat- 
terns also in position 4 hints at the presence of multiple 
similar channels throughout the inner wind. 

In C2, n e peaks at ~3000 cm -3 , in the south-west of 
position 1 as it overlaps with position 3 (i.e. the core of 
region A). This is distinctly different to what we see in 
CI: the western half of region A is much denser in C2, 
whereas the eastern half is considerably less dense. This 
type of result immediately highlights the advantage of 
spatia lly resolved observations over long-slit data (e.g. 
IW07d ). A secondary peak in C2 densities of just over 
1000 cm~ 3 is seen in the centre of region C (position 2). 
Again the morphology of C2 densities in region C differs 
significantly to those in CI. The C2 densities in the wind 
positions (4 and 5) follow the same pattern as seen in CI 
(falling off along the minor axis), but remain much higher 
for the same radii. The gas emitting C2 must therefore 
remain denser further out into the wind. It is worth 
noting here that the linear minor axis-oriented patterns 
are still visible, but appear to anti-correlate with those 
seen in CI. This is discussed further Section [5731 

The C3 density maps, tracing the density of the fainter 
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component of any split lines, are noisy and very incom- 
plete. There are, however, a number of points worth 
drawing attention to. C3 densities peak in the core of 
region C (at >2000 cm -3 , but with large uncertainties), 
approximately coincident with the flux peak. The gas 
knot in the centre of position 3 shows up clearly with 
densities between a few 100 and 1000 cm -3 . And in 
positions 4 and 5, C3 is detected in a number of north- 
west-south-east strips with densities of a few 100 cm~ 3 . 
This distribution, where the densities of the front- and 
back-sides of the expanding structures are clearly dif- 
ferent, emphasises that the structure of M82 is highly 
complex. 

To summarise this section, the CI densities do not fol- 
low the Ha emission morphology, but peak in the north 
of our spatial coverage, to the north of complexes A and 
C. This is coincident with the location of the major axis 
dust lane separatin g comple xes A and C from complexes 
D and E (see fig. 1 . IPaper H ). The C2 densities are much 
more consistent with the Ha flux distribution (i.e. the 
C2 densities peak within the SB complexes) , but remain 
much higher further out into the wind compared to CI. 
Since C3 represents the weaker component of any split 
lines detected, it is interesting that its density distribu- 
tion does not closely follow that of CI. Uncertainties on 
the C3 densities are, however, much higher than those 
for CI or C2. 

3.3. Excitation diagnostics 

The forbidden/rccomination line flux ratio of 
[Sn](A6717+A6731)/Ha can be used as an indica- 
tor of the number of ionizations per unit volume, and 
thus t he ionization p aram eter, U (|Veilleux fe Osterbrockl 
Il987t iDopita et afll2000L 12006( 1 . U is nearly constant 
within a single Stromgren sphere, a property that 
appears to hold over large scales in M82. [Sll]/Ha is 
also particularly sensitive to shock ionization because 
relatively cool, high-density regions form behind shock 
fronts which emit strong ly in [S ill thus producing a n 
enhancement in [Sn]/Ha l|Dopitalll997l : IOev et al.ll2000H . 
Maps of this ratio in the three line components are 
shown in Fig. [5] (left panels), where a single contour line 
at log[Sn]/Ha = —0.5 represents a fiducial threshold 
above which non-photoio nized emission is thought 
to play a dominant role (IDopita fc Sutherland! 119951 : 
iKewlev et all 120011: ICalzetti et al.l I2004H . In environ- 
mcnts such as this, the most likely excitation mechanism 
after photoionization is that of shocks. 

The ratio of [Nn]A6583/Ha can also be used 
in a similar fashion to [Sn]/Ha, as a diagnos- 
tic of the gas excitatio n level and to search for 
the presence of shocks (IVeilleux fc Osterbrockl 119871 : 
IDopita fc Sutherlandlll995H . In Fig. [51 (right panels), we 
show maps of this ratio for all three components; here 
the contour level represents a fiducial non-photoionized 
threshold of log[Nn]/Ha = —0.1. These maps are quali- 
tatively very similar to those of [S n]/Ha, and comparing 
the two figures gives an idea of the significance of the vari- 
ations and the uncertainties. To aid with this further, we 
have plotted the [Sn]/Ha and [Nil] /Ha points for each 
IFU position together in Fig. [101 separated into three 
graphs by line component. The dashed lines represent 
the aforementioned non-photoionization thresholds. 

The gas excitation in CI (Fig. [9l top panels) 



remains very uniform throughout the majority of 
positions 1, 2 and 3 (regions A and C), with 
mean values of log([S n]/Ha) = -0.96 ± 0.22, and 
log([Nn]A6583/Ha) = -0.28 ± 0.23. The errors quoted 
are the standard deviations of the samples, and their sim- 
ilarity shows just how constant these two ratios remain 
over these two positions. These results, in terms of both 
the absolute values and con stancy, are in good agreement 
with previous studies (IShopbell fc Bland-Hawthornl 
19981: iForster Schreiber et al.l 120011 : ISmith et alj 120061 : 
W07d). 

The C2 [Sll]/Ha and [Nil] /Ha line ratios arc both 
consistently lower than the average in two regions: com- 
plex A and the western half of complex C. For the most 
part these low excitation regions arc coincident with the 
corresponding C2 density peaks (Fig. [7|), suggesting that 
local changes in the ionization parameter play a strong 
part in setting the line ratios. 

In components CI and C2, high (non-photoionized) 
line ratios (both in [Sn]/Ha and [Nil] /Ha) are predom- 
inantly found in positions 4 and 5. This may indicate 
that, as the distance from the starburst complexes in- 
creases, shock excitation begins to play an increasing 
role in ionizing the gas. The phenomenon of increas- 
ing line ratios with radius has been found in many star- 
bursting galaxies on both the small and global scales 
dHeckman et al.lll995t lMartinlll997l : ICalzetti et aTlll999l . 
12004(1 . Where the line ratios are within the photoionized 
regime, the radial gradient can be explained as a dilution 
of the radiation field from the centralised source, but in 
many cases shock excitation is required to explain the 
ratios observed that are outside the reach of pure pho- 
toionization. 

Notably, one of the regions of high line ratios seen in 
the GMOS data, identifiable in both the C2 ma ps, is co- 
incident with the 'position 4 knot' identified in iPaperll . 
and discussed further in Section IBT41 The CI [Nil] /Ha 
map also exhibits very low ratios (log([N n]/Ha) < —0.6) 
in position 4 that are not coincident with any obvious fea- 
tures in the corresponding density map (Fig. [6]). In C3, 
high line ratios are found to the centre-east of position 
3 and to the north and east of the region C flux peak in 
position 2. 

In summary, the unprecedented spatial-resolution of 
these line ratio maps shows that (1) the excitation of 
the gas within the starburst complexes is extremely uni- 
form (and consistent with pure photoionization), despite 
the dramatic density gradients; and (2) regions of both 
high- and low-excitation gas are found within the inner 
wind flow, but only in very localised, distinct pockets. In 
CI, these regions do not correspond to any thing see n in 
any of the other maps presented here or in iPaper J . but 
the low-excitation regions in C2 are coincident with the 
density peaks. 

4. NEBULAR PROPERTIES OF THE DISK 

In this section we present and describe the 2D maps 
of Ha line flux and derived nebular properties for the 
DensePak dataset. In contrast to the small-scale, high- 
resolution GMOS-IFU observations of the nuclear re- 
gions, these data cover the majority of the whole cen- 
tral disk of M82 at conside rably low er spatial resolution. 
However, as highlighted in IPaper 1 the spectral resolu- 
tion and S/N are comparable, thus allowing us to make 
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TABLE 1 

Mean Ha component flux ratios for the GMOS and 
DensePak datasets. 



Ratio 


GMOS average 


DensePak average 


C2/C1 


1.65 1 


0.99 


C3/C1 


0.41 


0.43 



1 average excluding ratios >10 since these are likely to be spu- 
rious (see text). 

meaningful comparisons between the small- and largc- 
scales as sampled by the two datasets. 

4.1. Ha maps 

Fig. [11] shows the distribution of Ha flux in compo- 
nents 1, 2 and 3, and the equivalent width (EW; in A) 
of component 4, together with contours representing the 
WIYN i?-band image presented in Paper I. As expected, 
the flux map of CI (the bright narrow component) closely 
follows what is found through Ha narrow-band imaging: 
the bright nuclear starburst complexes stand out in posi- 
tion 1, whereas region B to the east of the disk (position 
3) is conspicuously faint in Ha. The broad component. 
C2, is detected over parts of all the positions, but is 
brightest in regions A and C and only detected in a few 
spaxels in region B. Over all the regions where C2 is de- 
tected it is, on average, of equal brightness to CI (see 
Table [l}. However, the C2/C1 flux ratio is not constant 
across the face of the disk, as can be seen in Fig. [121 
where we plot a map of the C2/C1 flux ratio. With this 
logarithmic colour stretch, it is clear that the ratio be- 
gins at low values near the disk midplane (0.3-0.5) and 
increases with height above and below the disk (up to 
-10). 

Component 3 (C3), the fainter narrow component, is 
detected over parts of the whole central starburst area, 
extending out to the north of the disk (position 4) and 
region B (position 3). As in the GMOS data, its intensity 
is on average —40% that of CI (TablcQJ. C4, the absorp- 
tion component, is observed over the whole of region B 
(in position 4 and a part of position 3) and to the south- 
west of cluster F in position 2. An absorption component 
is expected in these regio ns since they are l ocated out- 
side the current starburst (|Smith et a.l.ll2007h and hence 
contai n an older stellar population dominated by A-type 
stars (Ide Griis et al.ll20(n]jForster Schreiber et afll2003l : 
iKonstantopoulos et al.ll2008h . The Ha flux absorbed is 
at a maximum in region B, where EWs of —3 4 A are 
seen. An example line profile from region B is shown in 
fig. 5d of Paper I, clearly illustrating a broad and well- 
fitted absorption component. 

4.2. Electron densities 

Maps of the electron densities, rt c , measured from the 
[Sn] ratio (assuming T e = 10 4 ) are shown in Fig. [131 
Note, the velocities and FWHMs of the three [S n] com- 
ponents are consistent with those of Ha (Paper I) mean- 
ing that the components identified in both lines physi- 
cally correspond to one another, despite the lower S/N 
of the [Sn] line. In CI the highest densities are found 
in the eastern half of region A, peaking at —1000 cm -3 . 
The location of this peak agrees well with our GMOS 
data; its somewhat lower value can be attributed to the 
difference in spatial resolutions between the two datasets 



(i.e. the much finer spatial resolution of GMOS allows us 
to pick up the small-scale, localised regions of higher den- 
sity gas). The peak density in regi on C is lower (at 500- 
700 cm -3 ; in good agreement with lW0"7ch . and densities 
remain at a level of a few 100 throughout the central 
starburst region, particularly on the northern side. An- 
other peak is found in region B, with densities reaching 
500-700 cm- 3 . 

In C2 densities rise to >1000 cm -3 in the spaxels cov- 
ering the western half of region A and region E, in ex- 
cellent locational agreement with our GMOS data (again 
the lower absolute values can be attributed to the coarser 
spatial sampling of DensePak). Densities remain higher 
to much larger radii in C2, both in the major and minor 
axis directions. In order to highlight this further, we have 
plotted the radial distribution (relative to the nucleus) of 
both CI and C2 electron densities in the lower-right-hand 
panel of Fig. [T3j This clearly shows that the C2 densi- 
ties remain much higher out to larger radii relative to 
CI. This intriguing finding extends the trend we found 
in the GMOS data to the entire disk, and is discussed in 
r elation to the origin of C2 in S ection [5] 

lO'Connell fc Manganol (fl978h were the first to recog- 
nise that ionized gas densities in the starburst core of 
M82 were unusually high, finding [S ii]-derived electron 
densities of —1800 cm -3 in regions A and C, and sub- 
sequent multi-angle long-slit observatio ns showed that 
the g as densities rapidly fall off radially (jHeckman et al.1 
Il990t l . T he minor axis decrea se was accurately mea- 
sured by iMcKeith et alJ (|l995l ) using deeper long-slit 
data. They found electron densities of >1000 cm~ 3 near 
the nucleus (where they only detect their strongest line 
component), and by ±500 pc the densities had fallen 
to -100-200 cm" 3 . IW07c] presented HST/STIS spec- 
troscopy of the starburst core and the first evidence 
for how rapidly the densities vary with position: over 
only 2.5 arcsecs (45 pc) they found the de nsities to fall 
from —1800 cm -3 (near cluster Al; see also lSmith et al.l 
l2006t ) to -1000 cm" 3 in the centre of region A and 600- 
800 cm -3 in region C. 

Our DensePak measurements thus agree very well with 
these studies. The coarse spatial resolution of DensePak 
explains why we do not see the highest densities men- 
tioned in the literature and seen with our GMOS data. 
Having said that, our DensePak observations have the 
advantage of both of both wider two dimensional spatial 
coverage and high S/N, meaning we have been able track 
the density of both CI and C2 over a significant propor- 
tion of the inner disk. Finally, overall the densities of 
both CI and C2 stay higher on the northern side of the 
disk compared to that of the southern side. This can 
be understood as being a result of the 80° inclination of 
the disk, such that th e nort hern side is nearest to the 
observer (see fig. 13 of lW07d\ 

The main finding of this section is that although the 
the CI and C2 densities peak at approximately the same 
value in the centre of the starburst, the C2 densities re- 
main higher to much larger projected radii (i.e. both in 
the major and minor axis directions) than their CI equiv- 
alents. 

4.3. Excitation diagnostics 

Fig. [Hshows maps of the [S II] (A6717+A6731)/Ha and 
[N n]A6583/Ha line ratios. Since not all line components 
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were detected in each line, we have summed the flux of 
all identified components for each line before calculating 
the line ratios. 

A significant range in values are seen in both maps, 
with elevated ratios in region B (position 3 and the north 
east of position 4) and to the south-west of cluster F (po- 
sition 2). However, a closer inspection of the spectra from 
these regions reveals that the majority of spaxels with 
these high ratios exhibit an Ha absorption component 
- we have marked spaxels where an absorption compo- 
nent is detected with a cross in order to highlight the 
correspondence. Since absorption is only present in re- 
combination lines, the forbidden/recombination line ra- 
tios become biassed towards artificially high values when 
the absorption is present. This highlights the need to 
be careful of absorption effects in line ratios, particularly 
for ratios derived from images, or simple single Gaussian 
fitting, where absorption cannot be accounted for. 

Disregarding these Ha-absorption affected regions, 
the two line ratios remain fairly uniform (particularly 
[Nil] /Ha). The mean ratio of log([Nn]/Ha) for all spax- 
els not containing absorption is —0.3 ±0.1, and similarly 
for log([S n]/Ha) is —0.8 ± 0.3. These values are in very 
g ood agreement with the GMOS data (Section R?~3")) . 

10'Connell fc Manganol (fl97l find a log([N Il]/Ha) 
ratio for the nuclear regions (including com- 
plexes A, C and E) of — -28 f rom long-slit data. 
IShopbell fc Bland-Hawthorn! (|1998D find uniform and 
low values of log([N n]/Ha) in the inner 1 kpc of the 
M82 outflow ranging from —0.5 to —0.2. Both these 
results are consistent wit h our measurements, and in 
very good agreement with IW07d. T he line ratio maps of 
IShopbell fc Bland-Hawthornl (|1998fl also show very high 
values further out in the disk (to the east and west of the 
nucleus), in the same place as we find exaggerated values 
due to Ha abs orption. At minor axis height s of 1 k pc or 
more (>60"), IShopbell fc Bland-Hawthornl (fl998h find 
increasing values of [Nil] /Ha which they attribute to a 
dilution of the radiation field in the outer-wind together 
with a possible increasing influence of shocks. However, 
our data do not reach this far in radial extent so we do 
not see this effect. 

5. DISCUSSION 

Using all the information we now have about the gas 
conditions within the M82 starburst, we can begin to put 
together a model of the ISM in the starburst core, and 
examine how these conditions change as we move out of 
the disk into the wind. 

5.1. ISM conditions within the starburst core 

The electron densities we have mapped peak in the 
starburst core, fall off dramatically as we leave the disk, 
and exhibit significant local variations. The peak val- 
ues (few 1000 cm -3 ) imply very large interstellar pres- 
sures (P/k > 10 7 ). Although these results have been 
found before by a number of independent studies (e.g. 
iHeckman et all HTfflOl : iMcKeith et all [1991 IW07d\ our 
high-resolution, two-dimensional maps show that the 
peaks are very localised and that there is a large amount 
of structure within the density distribution in the star- 
burst core. 

Contrary to expectation, the CI densities peak not in 
the starburst complexes themselves, but to the north, in 



the major axis dus t lane. Recalling one of our conclu- 
sions from iPaper II - that the emission of CI is weighted 
towards the inner-most regions of the starburst (i.e. the 
base of the outflow) - we see that these density results 
support this idea, and indicate that the body of bright, 
relatively quiescent gas from which CI is emitted, is not 
associated directly with the main starburst complexes 
themselves, but extends well beyond the complexes into 
the dense dust lane. C2, however, is at its densest within 
complex A (and to a lesser extent C); its overall distribu- 
tion follows the flux morphology much more closely than 
CI. Its peaks are also of a few 1000 cm -3 implying large 
pressures. Th ese resul ts are again in support of our con- 
clusions from IPaper 1 where we suggested that the C2 
emission is weighted towards regions further out along 
the line-of-sight. 

In this paper we have also presented evidence at high 
spatial resolution for the impressive uniformity of the line 
ratios and excitations within complexes A and C (par- 
ticularl y in component CI), firs t high light ed by studies 
such a s IForster Schreiber et al.l (|2001[ ) and iBeirao et al.l 
(|2008h . We have also confirmed how this uniformity can 
be extended to the whole inner disk region covered by 
the DensePak fields-of-view. The unprecedented spatial- 
resolution and spatial coverage of the line ratio maps 
shows that (1) the excitation of the gas within the star- 
burst complexes is extremely uniform (and consistent 
with pure photoionization) , despite the dramatic den- 
sity gradients; and (2) regions of both high- and low- 
excitation gas are found within the inner wind flow, but 
only in very localised, distinct pockets. 

We will now use these density and excitation re- 
sults, together with work from the literature, to build 
up a model to describe the ISM conditions within 
the starburst core. As described in the introduc- 
tion, the modelling of observations of gas in the cen- 
tral regions made from the opticalradio indicate a frag- 
mented (shredded/fractious) ISM, consisting of a range 
of clouds from small/dense clumps with low filling factors 
(<1 pc, n e > 10 4 cm" 3 ) to larger filling factor, less dense 
gas (e.g. lLord et al.lfl996l : IForster Schreiber et al1l2001l : 
iFuente et al.ll2008ft . These are bathed with an intense ra- 
diation field and embedded in an extensive high tempera- 
ture (T > 10 6 K), X -ray-emitting ISM. Emission measure 
results from lW07c1 suggest that the most compact clouds 
are found in the clump cores where the star formation 
is most intense and the gas pressures are high, whereas 
the cloud sizes in the inter-complex regions are larger. 
The observed near-constant state of the starburst (ion- 
ization parameter and line diagnostics) must therefore 
be a consequence of the small cloud sizes, allowing the 
gas conditions to respond quickly to any changes. This 
is rather akin to the thcrmalisation of a large number 
of individual particles in a plasma. To this picture, we 
are now able to add that the denser clouds are not only 
bathed in UV but are also interacting with an extensive, 
very hot, and presumably turbulent ISM phase that is 
responsible for the X-ray emission. The left-centre panel 
of Fig. [TJ5] contains a schematic representation illustrat- 
ing these ideas. Since we associate the broad component 
C2 with emission from turbulent mixing layers on the 
surface of interstellar clouds, the existence of so many 
small, compact, high density clumps embedded in a hot, 
fast-flowing wind must make for an ideal environment 
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for the copious emission of broad lines. These surface 
wind-clump interaction layers must also be where mass 
is evaporated and stripped off the clouds by the intense 
radiation field and wind flows. This suggests that the 
ma jority of the mass loaded into the wind must origi- 
nate in these clumps, where the density of the gas being 
stripped is highest. This argues strongly against C2 orig- 
inating from scattered light from the nucleus. A further 
implication of the existence of so many dense clouds in 
complexes A and C is that there must still be a significant 
gas supply within the starburst core, despite the intense 
star formation, heavy mass-loading and wind outflows. 

In our model, C2 is only emitted from turbulent mix- 
ing layers on cloud surfaces in direct line-of-sight with a 
wind source, whereas CI can be emitted from any 10 K 
gas located within the Stromgrcn radius of the starburst. 
This includes the static ionized layers beneath the TMLs 
(bottom panel of Fig. [T6|) , ionized gas clouds which are 
not in a direct line-of-sight with an ionizing/wind source 
(i.e. wind shadowed - ionizing radiation is easily scat- 
tered an example of a cloud in a wind shadow is illus- 
trated in the bottom-left corner of the centre-left panel 
in Fig.fTH]). and general diffuse ionized gas. Furthermore, 
the excitation of the gas within the starburst complexes 
is extremely uniform, and consistent with pure photoion- 
ization. This can be explained through the well-mixed 
nature of the clusters and clouds, which transforms the 
presumably very inhomogeneous distribution of ionizing 
photons from the clusters into a uniform escaping radi- 
ation field. A helpful analogy might be to liken the sur- 
faces of the starburst complexes to the surface of the Sun, 
where there are many gas phases (temperatures and den- 
sities), a strong density gradient, chaotic dynamics and 
morphology, and denser pockets of gas in the form of 
loops and filaments (shredded shells/bubbles) entrained 
in the outflow/wind, all being accelerated outward. 

As a final note to this section, we turn our atten- 
tion to the the ~7 Myr old super star cluster M82-A1 
([Smith et alj|2006t) . In that study we found M82-A1 to 
have followed a non-standard evolutionary path for a star 
cluster, since it is surrounded by a lar ge, pressurised H n 
nebula that is stifling its wind (see also lSilich et al.|[2007f ) . 
If the H II shell forms a full sphere around the cluster and 
this stage is not simply transient, then it is difficult to 
explain how the superwind is powered if one assumes 
M82-A1 to be representative of a typical cluster in re- 
gion A. However, our GMOS maps show that M82-A1 is 
in fact embedded in one of the highest CI density peaks 
in the starburst core. Recalling that the inclination of 
the disk is such that on the southern side we are seeing 
up into the starburst (the northern side is par tially ob- 
scured by the front of the disk; see e.g. IW07c[ ). we can 
now see that the rapid fall in densities from the north of 
our IFU fields into the cores of complexes A and C shows 
that the conditions towards the south of the complexes 
are much more favourable to the escape of the individual 
cluster winds. This implies that clusters towards the out- 
side surfaces of the complexes may in fact be the primary 
drivers of the superwind. 

5.2. Inner wind 

As we move away from the plane of the disk, there 
are a number of results that indicate that emission from 
C2 becomes increasingly more important (compared to 



CI). Firstly the C2/C1 flux ratio increases (all the 
way through our GMOS fields out to the extent of 
our DensePak data; Sections 13.11 and 14. ip . Since the 
C2/C1 ratio reflects the proportion of the volume occu- 
pied by TML vs. standard photoionized gas, this indi- 
cates that within the energy injection zone (thought to 
be defined by the minor axis velocity inflection points ; 
McKeith et atl[T995l : iShopbell fc Bland-HawthonJ 11991 
W07d) the turbulent motions must play an increasing 
role with radius. Secondly, the electron densities in C2 
remain higher to larger radii (Sections 13.21 and | 4.2[) . sug- 
gesting that away from the dense core of the galaxy, then 
density in the TMLs is preferentially higher than in the 
narrow line-emitting gas. Although the majority of the 
mass loading might take place within the starburst com- 
plexes, the region of mass loading must extend far out 
into the inner wind. Thirdly, the line widths of C2, al- 
though already high in the nuclear r egions, also increase 
with height a bove the disk ([Paper J) , peaking at the lo- 
cation of the iMcKeith et al.l (j!995l )' minor axis velocity 
inflection points. Finally, it is possible that these effects 
are also partly due to the weakening or dilution of the 
CI emission with height. What happens to these rela- 
tionships beyond the energy injection zone (beyond the 
minor axis velocity inflection points) is something that 
should be investigated. 

On average the C2/C1 flux ratios are consider- 
ably higher than what w as found in NGC 1569 
([Westmoquette et al.ll2007bT ). Since C2 is produced at 
the interface layers of interstellar clouds, the significantly 
larger gas mass and number of clusters in M82 has re- 
sulted in C2 being more prominent and may mean that 
a larger percentage of mass is being loaded into the out- 
flow. 

5.3. Outflow channel 

In iPaper II we identified a discrete outflow channel 
in IFU position 5 characterised by expanding, narrow- 
line-cmitting walls surrounding a cavity of broad-linc- 
emitting gas. With the additional results that we have 
presented here, we can explore the characteristics and 
implications of this outflow channel further. 

Fig.[15]shows Ha position- velocity plots extracted from 
seven evenly spaced (every 2") pseudo-slits positioned 
parallel to the major axis, stepping down through GMOS 
IFU position 5. Also plotted are the CI and C2 electron 
densities extracted from the same pseudo-slits (the C3 
densities have been omitted for clarity, but are consis- 
tent with those of CI; see Figs. [5] and [5]). A coherent, 
expanding cylinder, as indicated by the overlaid dashed 
ellipses, can be traced in the narrow-line (CI and C3) 
kinematics extending along the whole length of position 
5. These velocity ellipses show that the structure, while 
>7" (120 pc) in length, increases slightly in physical 
width from the north to south from ~2-3" (35-50 pc), 
and maintains an approximately constant (subsonic) ex- 
pansion velocity of ^60 kms -1 . This is indicated in a 
schematic cartoon of the outflow channel in the right- 
hand panel of Fig. [16] The velocities of the broad C2 
line suggest that it originates within the space between 
the two walls, with central velocities at rest with respect 
to the C1-C3 expansion. The [S ii]-derived electron den- 
sity plots demonstrate how the walls of the channel are 
defined by peaks in the C2 densities (700-1000 enr 3 ) 
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and troughs in the CI densities (~200 cm -3 ), whereas 
the centre of the channel is defined by higher CI den- 
sities (>500 cm -3 ) that decrease with radius. What is 
not shown in this plot is that the widths of C2 appear 
to be broader (~250 kms -1 ) near the channel walls and 
narrower (>150 kms -1 ) in the centre (jPaper J ). 

Thus, in our model the channel cavity contains the 
hot outflowing wind fluid. As it flows outwards it in- 
teracts with the cooler, denser walls of the channel and 
with clumps of entrained material to produce broad-line 
emission, while the walls themselves emit primarily the 
narrow lines CI and C3. Finally, we note that the lin- 
ear minor axis-oriented patterns observed in th e density , 
FWHM and radial velocity maps (Section 13.21 IPaper J ) 
hint at the existence of multiple similar channels within 
the sout hern outflow, as is predict ed by models such as 
those of ITenorio-Tagle et all (|2003t) . 

The presence of organized channels which funnel shock- 
heated gas outwards from starburst complexes has a 
number of important implications. These structures re- 
semble scale d-up versions of the cluste r wind channels 
described by ITenorio-Tagle et al l (|2003t) . and as in their 
models, will tend to collimate hot gas outflows. The col- 
limated outflows are a possible source for the well defined 
loop structures of ionized gas seen in Ha emission above 
the main starburst complexes in M82. Indeed hints of 
the cha nnel can be seen in the Ha morphology (see e.g. 
fig 1 in IPaper If ) in the form of enhanced emission fol- 
lowing the channel walls. Without the line ratio and 
kincmatical evidence, however, it would have been im- 
possible to distinguish from any other formation in this 
chaotic region. Channelled flows will also influence the 
mass loading process, and it seems likely that some of the 
gas within such columns could experience minimal mass 
loading. In this case they will improve the efficiency in 
launching a galactic wind in that some hot gas can flow 
out of the disk relatively unimpeded. In addition out- 
flow channels can provide low opacity pathways for the 
escape of UV radiation from starburst complexes which 
in turn can help ionize the large quantities of extraplanar 
gas associated with this galaxy. 

5.4. Position 4 knot 

In Paper I we identified a small region in CMOS IFU 
position 4, coincident with a bright, point-like object 
on HST/ ACS imaging, with distinctly blucshiftcd (by 
~140 kms" 1 ) and broad (<350kms -1 ) Ha emission that 
we referred to as the "position 4 knot" . From the data 
presented in this paper, we see that this knot does not 
show up in either the flux or density maps. However, 
a small number of spaxels near the location of the knot 
exhibit a slight enhancement in [Sii]/Ha C2 and a sig- 
nificant enhancement in [Nil] /Ha C2 compared to the 
surroundings (Fig [9]). To confirm these results, we ex- 
tracted and summed the spectra from 10 spaxels centred 
on the knot; the resulting spectra and multi-Gaussian 
line fits are shown in Fig. [T7l The broad (FWHM ~ 
330 kms" 1 ), blucshiftcd (—140 kms -1 ) component is 
clearly identifiable in all five emission lines. The line 
component properties are given in Table [5J The elec- 
tron densities derived from the [S n] line ratios indicate 
that the gas from which the broad component is being 
emitted is denser than the narrow line gas (~1000 cm -3 
compared to ~500-900 cm -3 ). The line ratio enhance- 



ments described above do not show up in the summed 
spectrum since they have been diluted by the surround- 
ing spaxels. 

Inspection of archived HST / ACS HRC images (Prop 
ID: 10609, PI: Vacca) reveals one other similar nearby 
point-like source (hereafter source 2), the location of 
which, together with the position 4 knot, is identified on 
a F330W (~J7), F435W (~B), F550M (~V) and F814W 
(~/) HRC colour composi te in Fig. IT51 Measurement of 
their sizes (using ishape; lLarsenir2004l ) show that they 
are both unresolved (FWHM < 0'.'06) in each of the 
ACS/HRC images, including F658N (-Ha). For both 
objects, we find no evidence for variability (>0.2 mag) 
in any of the aforementioned filters over the epoch range 
September 1995 to March 2006 (the range covered by 
all the HST WFPC2 and ACS images available in the 
archive). Absolute F555W magnitudes (My) and colours 
are given in Table [3l No evidence of the object is found 
in any of the available ffST/NICMOS NIC1/2 images 
(F160W, F164N, F187N, F190N). 

Considering our spectroscopic and photometric evi- 
dence and their unresolved nature in the HST imaging, 
there are a number of plausible explanations for these 
sources. We can immediately dismiss the position 4 knot 
as being a background object based on its radial velocity, 
it only having a difference of ~140 kms -1 relative to the 
systemic velocity of M82. Furthermore, since it does not 
exhibit an enhanced density compared to its surround- 
ings (Fig. [7]), does not show up in any of our IFU line 
or continuum flux maps, and does no t correlate with any 
known radio source at this location dHuang et al.l 119941 : 
IWills et al.|[i~997l ; iMcDonald et al.ll2002f) . it is unlikely to 
be a supernova remnant. 

We find the most likely explanation is that the position 
4 knot is a hot stellar object immersed in an unresolved 
ionized circumstcllar nebula. The FWHM of the broad 
component then represents the expansion of the nebula 
of —165 kms -1 . We derive an absolute magnitude of 
My ~ —5.2 which increases to ——7 if we assume a V- 
band extinction of ~2 mags, suggesting a massive star. 
The photometry shows that the star is fairly red but 
this may be due to a dusty circumstcllar nebula. The 
high [N n] /Ha ratio suggests that the nebula may be en- 
riched in nitrogen. All these findings lead us to speculate 
that the position 4 knot may be a luminous b lue vari- 
able (LBV) star (jHumphrevs fc Davidson! I1994T ). These 
objects are usually surrounded by nitrogen-rich circum- 
stellar material which has been ejected as a result of an 
unstable event in the past. These stars can spend many 
years in a hot, quiescent phase and have similar abso- 
lute magnitudes to the position 4 knot . The large radial 
velocity of the emission from this object is intriguing; it 
might indicate that the star was ejected from the cluster 
in which it was formed, thus explaining how such a star is 
found so far from the main star forming complexes. Mea- 
suring a projected distance of the star from the nearest 
star- forming region is ~5-6" (~90-100 pc), and even as- 
suming a conservative estimate of the age of the LBV (a 
60 Mq mass LBV would be <4 Myr old), the predicted 
velocity is quite plausible (<50kms -1 ). 

Source 2, for which we only have photometry mea- 
surements, is brighter and redder than the position 4 
knot (Table [3]). Although it is also unresolved on the 
HST / ACS HRC images and also does not correlate with 
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TABLE 2 

Emission line Gaussian component measurements from the summed spectrum (Fig. [TTIi of the position 4 knot. 



Component 


Ha FWHM 


Ha radial velocity 


log([N II] /Ha) 


log([Sll]/Ha) 


Density 




(kms -1 ) 


(kms- 1 ) 






(cm -3 ) 


CI 


156 


40 


-0.34 


-0.59 


480 


C2 


333 


-138 


-0.43 


-0.64 


1100 


C3 


53 


37 


-0.42 


-0.81 


900 



TABLE 3 

Point source photometry, where U = F330W, B = F435W, V = F555W, I = F814W. 



Object 




Coordinates 




V 


M v 


U-V 


B — V 


V -I 


B — I 


U -I 


U-B 






(J2000) 










— (mags) — 








1 (pos4 knot) 


09 h 55 m 


54!06 69° 40' 


44'.' 35 


22.16 


-5.22 


1.14 


0.98 


1.40 


2.38 


2.45 


0.16 


2 


09 h 55 m 


53?03 69° 40' 


39'.' 60 


20.46 


-6.92 


0.68 


0.62 


0.79 


1.41 


1.47 


0.06 



any known radio sources. Spectroscopic measurements, 
well within the capabilities of 8-10 m class telescopes, 
are needed to ascertain the true nature of these objects. 

6. SUMMARY 

By combining both detailed (~0'.'7) observations of 
M82's starburst core and inner wind made with the Gem- 
ini GMOS-IFU, and wider-field observations of the inner 
~2 kpc of the disk made with the WIYN DensePak in- 
strument, we have been able to build up a comprehensive 
picture of the state of the ionized gas within this galaxy. 
To do this we split the study, first focussing on the ion- 
ized gaseous and st ellar dyn amics of the M82 starburst 
core and inner disk (|Paper J ), then on the derived nebu- 
lar properties of the ionized gas (this paper). 

The dynamics showed that the ionized gas within the 
starburst is dynamically highly complex with much lo- 
calised line splitting, and many overlapping expanding 
structures located at different radii. An offset of ^12° 
between the rotation axis of the stars and gas provides 
further evidence of the effects of the interaction with M81 
that is thought to have initially triggered the starburst. 
We have identified an underlying broad component in the 
optical nebular emission lines throu ghout the inner disk, 
and fol lowing the re cent work by IWestmoquette et"aLl 
(|2007al lbl) and IW07d . have associated it with turbulent 
layers on the surfaces of interstellar clouds set up by the 
impact of high energy photons and fast-flowing winds 
from the embedded young star clusters, and hence the 
sites of mass loading within the wind outflow. 

Maps of the [S ii] -derived electron densities show that 
in both the narrow and broad component gas, signifi- 
cant small-scale local variations are present, and that 
the peaks are of a few 1000 cm -3 implying very large in- 
terstellar pressures (P/k> 10 7 ). Despite these variations 
and a dramatic fall-off in thermal pressure with height, 
we still find an impressive uniformity in the line ratios 
and excitations within the whole inner disk region, mean- 
ing that the excitation of the gas within the starburst is 
extremely uniform. That the ratios are consistent with 
pure photoionization suggests that shocks must not play 
an important role in setting the gas ionization within the 
disk. Our results thus join with those of previous studies 
in showing that the inner region of M82 has many of the 
characteristics of a single supergiant Stromgren sphere. 

These findings, when put together with the body of lit- 



erature built up over the last decade on the state of the 
ISM in M82, imply that the ISM is highly fragmented 
into a range of clouds from small/dense clumps with low 
filling factors (<1 pc, n c > 10 4 cm -3 ) to larger filling fac- 
tor, less dense gas. The most compact clouds seem to be 
found in the clump cores, whereas the cloud sizes in the 
inter-complex region are larger. These dense clouds are 
bathed with an intense radiation field and embedded in 
an extensive high temperature (T > 10 6 K) ISM that is 
a product of the high star formation rates in the star- 
burst zones of M82. The near-constant ionization state 
of the ~10 4 K gas throughout the M82 starburst zone can 
be explained as a consequence of the small cloud sizes, 
which allow the gas conditions to respond quickly to any 
changes. 

In our model, schematically represented in the left- 
centre panel of Fig. [16j the broad component (C2) is 
only emitted from turbulent mixing layers on cloud sur- 
faces in the direct line-of-sight of a wind source. The 
large number of compact clouds mixed in with the wind 
sources of course provides an ideal environment for copi- 
ous emission of these broad lines. In contrast, the narrow 
components are emitted by relatively quiescent 10 4 K gas 
within the M82 Stromgren sphere, including photoioniza- 
tion fronts beneath the TMLs on the surfaces of dense 
clouds, ionized gas clouds which are not in the direct line- 
of-sight of an ionizing/wind source (i.e. wind shadowed), 
and general diffuse ionized gas. 

Previous obse rvations of the brig ht region A super star 
cluster M82-A1 ([Smith et al.ll2006h indicate that its wind 
has been stalled by the high interstellar pressures sur- 
rounding it. If this cluster is surrounded by a static 
spherical shell then it would be difficult to undertsand 
how it, and possibly other massive star clusters, con- 
tribute to launching the M82 supcrwind. However, we 
can now see that the rapid fall in densities from the 
north of our IFU fields into the cores of complexes A 
and C and beyond shows that the conditions towards 
the south of the complexes (their outward surfaces) may 
be much more favourable to the escape of the individ- 
ual cluster winds. Clusters towards the outside surfaces 
of the starburst complexes may in fact be the primary 
drivers of the supe rwind, not clusters such as M82-A1 
([Smith et all l2006t ) which our maps show is embedded 
in one of the densest parts of the starburst. This idea 
is consistent with the asymmetric distribution of extra- 
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planar material in M82: a preponderanc e of blue source s 
are found to the south of of the disk by iDavidgd (|2008f ) , 
and a clear ly dominant southe rn blowout region is seen 
in the UV (jHoopes et al.ll2005f) . 

By extr acting m ajor axis position-velocity plots, we 
showed in iPaper II that the broad-line emitting gas ro- 
tates, on average, at a slower velocity than the narrow- 
line gas. The identification of a discrete channel within 
the outflow allowed us to explain this velocity differ- 
ence: if the broad component, representing turbulent 
gas stirred up by the interaction between the hot, fast 
wind and cooler gas, becomes more dominant as the gas 
becomes increasingly perturbed and entrained along the 
length of the channel, then its emission will be weighted 
towards regions further out along the flow compared to 
the narrow component, whose emission is weighted to- 
wards the base of the outflow channel. 

This idea is supported by the findings of this paper, 
where we saw that the CI densities peak not in the star- 
burst complexes themselves, but further in within the 
major axis dust lane, whereas C2 is at its densest within 
complexes A and C (i.e. further out) and exhibits a slower 
fall off (with height) than CI. This result, when com- 
bined with the fact that the C2/C1 flux ratio increases, 
and the C2 line widths, although already high in the 
nuclear regions, also increase (peaking at the location 
of the minor axis velocity inflection points) imply that, 
within the energy injection zone (as defined by the minor 
axis velocity inflection points), the turbulent motions as 
traced by the broad C2 component appear to play an 
increasing role with height. This behavior parallels the 
situation that we found in NGC 1569, and is likely to 
result from the operation of similar physic al processes in 
both galaxies (jWestmoquette et al.l 12008'). If the broad 
component indeed originates from the wind-clump inter- 
action layers where material is also ablatively stripped 
from the clouds, then although the majority of the mass 
loading might take place within the starburst complexes, 
the region of mass loading must extend far out into the 
inner wind. 

The narrow-line (CI and C3) Ha position-velocity 
plots extracted from CMOS IFU position 5 have allowed 
us to look in more detail at the aforementioned outflow 
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channel. The channel appears as a coherent, expanding 
cylindrical structure extending along the whole length of 
position 5, increasing slightly in physical width from the 
north to south from ~2-3" (35-50 pc), and maintaining 
an approximately constant (subsonic) expansion velocity 
of ~60 kms -1 . Velocities of the broad C2 line suggest 
that it originates within the space between the two walls. 
The walls of the channel are defined by peaks in the C2 
densities (700-1000 cm -3 ) and troughs in the CI den- 
sities (~200 cm -3 ), whereas the centre of the channel 
is defined by higher CI densities (>500 cm -3 ) that de- 
crease with radius. In our model, the channel cavity 
contains the hot wind fluid that, as it flows outward, in- 
teracts with the cooler, denser walls of the channel and 
with clumps of entrained material to produce broad-line 
emission, while the walls themselves emit primarily the 
narrow lines CI and C3. Hints are seen that there may 
be multiple similar channels within the outflow. 

Finally, we have argued that a point-like knot 
identified in CMOS position 4, exhibiting blueshifted 
(~140 kms -1 ), broad (<350 kms" 1 ) Ha emission and 
enhanced [Sii]/Ha and [Nil] /Ha ratios, is most likely a 
luminous blue variable (LBV) type object, that may have 
been ejected from the cluster in which it was formed. 
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Fig. 1.— HST F658N (Ho) image of M82 showing the positions of the C MOS (white) and De nsePak (black) IFUs. Some of the 
bright starburst co mplexes /clusters are marked with letters (nomenclature from O'Conncll et al. 1995), the position of the 2.2 (im nucleus 
lILester et al.|[T990h is shown with a white cross, and the galaxy major axis (PA = 65°) is indicated by a solid line. 
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Fig. 2. — Ha flux maps for CI. A scale bar is given for each map in arbitrary (but relative) flux units. Contours of the ACS F435W 
(B-band) are overplotted to illustrate the locations and extents of regions A and C. The maximum size of the seeing disk (0"8) is indicated 
by the black circle in the lower-left and applies to all subsequent CMOS figures. 




arcsecs 

Fig. 3. — Ha flux maps for C2. Labels as for Fig. [2] 




arcsecs 

Fig. 4. — Ho flux maps for C3. Labels as for Fig. [2] 




arcsecs 



Fig. 5. — GMOS Ha C2/C1 flux ratio maps. The single contour level represents a ratio of 1 (i.e. C2 = CI), as also indicated with a 
marker in the scale bar. (A color version of this figure is available in the online journal.) 
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Fig. 6. — Maps of the [S ll]AA6717,6731 derived electron densities (n e , in units of cm -3 ), for CI. The lower limit to the colour range 
equates to 100 cm -3 , the low-density limit of this indicator, and the contours represent 1000, 2000, 3000, and 4000 cm -3 , as marked on 
the scale bars. (A color version of this figure is available in the online journal.) 
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Fig. 7. — Maps of the electron densities for C2. Labels as for Fig. [6] (A color version of this figure is available in the online journal.) 
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Fig. 8. — Maps of the electron densities for C3. Labels as for Fig. [6] (A color version of this figure is available in the online journal.) 
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Fig. 9.— Maps of the flux ratio of log([S ll](A6717+A6731)/Ho) (loft panels) and log([N II] A6583/Ha) (right panels). Contour lines are 
plotted at —0.5 ([Sll]/Ha maps, as indicated in the colour-bar) and —0.1 ([Nil] /Ha maps, also indicated in the colour-bar), representing 
fiducial ratios above which the excitation is likely to be dominated by non-photoionizing processes (see text). (A color version of this figure 
is available in the online journal.) 
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Fig. 10.— Plots of log([S n](A6717+A6731)/Ha) vs. log([N ll]A6583/Ha) for each IFU position, separated by line component. The vertical 
dashed lines represent the non-photoionziation thresholds from Fig. [9] (A color version of this figure is available in the online journal.) 
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Fig. 11. — Line flux maps for the four identified Ha components. The CI, C2 and C3 maps are in arbitrary (but relative) flux units, 
whereas the C4 (the absorption component) map is in units of equivalent width (EW) in A. Ovcrplotted on the CI map are contours from 
the log scaled WIYN ii-band image (see Paper I) to indicate the location of the starburst complexes. The origin is set to the position of 
the 2.2 (im nucleus. (A color version of this figure is available in the online journal.) 




Fig. 12. — DensePak Ho C2/C1 flux ratio map (in logarithmic units). As in Fig. [5] the marker in the scale bar indicates a ratio of 1 (i.e. 
C2 = CI). (A color version of this figure is available in the online journal.) 
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Fig. 13. — Maps of the [So] derived electron density, ra e , in the four DensePak fields. Units are in cm 3 ; markers in the scale bar indicate 
500, 750, 1000 and 1250 cm~ 3 . Uncertainties lie in the range ±500 for the highest densities to ±50 for the lowest. Ovcrplotted on the CI 
and C2 maps are contours from the log scaled WIYN ij-band image (see Paper I) to indicate the location of the starburst complexes. The 
lower-right-hand panel shows a radial plot of the electron densities of CI and C2, highlighting the more extended distribution of higher C2 
densities compared to CI. (A color version of this figure is available in the online journal.) 
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Fig. 14. — Flux ratio maps of log([S ll](A6717+A6731)/Ho) (upper panel) and log([N II] A6583/Ha) (lower panel) - note the different 
colour bar scales. The ratios were calculated after summing the flux from all identified components of the lines. Spaxels marked with crosses 
are those that contain an identified Ha absorption component, illustrating how the absorption is biassing the line ratios to artificially high 
values. (A color version of this figure is available in the online journal.) 
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Fig. 15. — Position-velocity plots for seven evenly spaced (every 2") pseudo-slits positioned parallel to the major axis, stepping down 
through GMOS IFU position 5 (i-iv = north-south). Ha CI is plotted with black circles, C2 with red squares, and C3 with green stars. 
The dashed ellipses highlight the expanding outflow channel structure traced by Ho CI and C3; the velocity of the broad Ho C2 indicates 
that it originates from within the channel. Also plotted are the corresponding CI (solid line) and C2 (dashed line) electron densities (C3 
has been excluded for clarity but is consistent with CI), showing the anti-correlation in CI and C2 densities within the channel (CI densest 
in the channel centre and falls off with radius, whereas C2 is densest in the channel walls). (A color version of this figure is available in 
the online journal.) 
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Fig. 16. — Cartoon model depicting the ideas developed in this work regarding the state of the ISM in M82. At the top is a representation 
of the disk showing the location of starburst complexes A and C, the nucleus (white cross), and the major axis dust lane (hatched region). 
The left-hand panel shows schematic model of the ISM in the core of one of the starburst complexes. The gas is fragmented into a large 
range of cloud sizes, the star clusters are well mixed in with the clouds, and both are embedded in a hot medium at very high pressure. 
The characteristic size-scales of the clouds and the clusters are similar (few pc). Turbulent mixing layers (TMLs) form on the surfaces 
of the clouds directly facing the oncoming winds from the clusters, whereas the surfaces of all the clouds are fully ionized by the energy 
supplied by the star clusters, even if they are located in "shadowed" regions (an example of a cloud in a wind shadow is illustrated at 
the bottom-left). The bottom panel shows a single cloud b e ing ir radiated by ionizing photons and impacted by the fast-flowing wind 
from a nearby star cluster. Following the idea of Bincttc ct al. (2009), a stratified cloud structure forms, with a TML on the cloud surface 
facing the cluster (from which the broad line is emitted), followed by a static photoionized layer (narrow lines). The core of the cloud is 
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Fig. 18.— HST ACS/HRC colour composite (F330W, F435W, F550M, F814W) image with the position 4 knot and the other similar 
objects identified, together with the GMOS IFU FoVs. (A color version of this figure is available in the online journal.) 



